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Earlier  results  for  coherent  propagation  of  light  in  correlated  random  distributions  of  dielectric  particles  of  radius 
a  (with  minimum  separation  b  >  2 a  small  compared  with  wavelength  X  =  2tr/6)  are  generalized  to  obtain  the  re¬ 
fractive  and  absorptive  terms  to  order  fka  )2.  The  present  results  include  the  earlier  multiple  scattering  by  electric 
dipoles  as  well  as  scattering  and  multipole  coupling  by  magnetic  dipoles  and  electric  quadrupoles.  The  correlation 
aspects  are  determined  by  the  statistical-mechanics  radial  distribution  function  /(R)  for  impenetrable  particles 
of  diameter  b.  The  new  terms  for  slab  scatterers  and  spheres  involve  the  integral  of  fR  (first  moment)  or  of  f  In  R 
for  cylinders.  The  new  packing  factor  is  evaluated  exactly  for  slabs  as  a  simple  algebraic  function  of  the  volume 
fraction  w.  and  it  is  shown  that  the  bulk  index  of  refraction  reduces  to  that  of  one  particle  in  the  limit  w  =  1.  A  sim¬ 
ilar  result  is  achieved  for  spheres  in  terms  of  the  Percus-Yevick  approximation  and  the  unrealizable  limit  w  =  1. 


INTRODUCTION 

Earlier  papers1-  developed  simple  forms  for  the  coherent  bulk 
index  of  refraction  ( rj  =  y  e)  for  correlated  random  distribu¬ 
tions  of  dielectric  particles  with  minimum  separation  ( b )  of 
centers  small  compared  with  wavelength  (X  =  k/2tr).  Writing 
the  index  as  r)  =  ty  +  ir\a  +  ir]s,we  applied  general  scattering 
theory1-4  to  the  range  of  small  kb  to  obtain  results  for  the  re¬ 
fractive  (rjr)  and  absorptive  (tj„)  terms  that  were  explicitly 
independent  of  X  and  to  obtain  corresponding  results  for  the 
scattering  (»>,)  loss  term  to  lowest  order  in  X.  The  explicit 
approximations2  for  rjr  and  rja  for  spheres,  cylinders,  and  slabs 
(m  =  3.  2.  1,  respectively)  depended  only  on  the  particles’ 
radius  or  half-width  (a),  their  complex  index  of  refraction  (rj' 
=  \  T7),  and  their  average  number  ip)  per  unit  volume;  they 
exhibited  the  statistical  aspect  of  the  problem  only  in  the 
volume  fraction  w  =  pv,  with  c  =  t'la)  as  the  volume  of  one 
particle.  The  corresponding  scattering  terms  17,  were  addi¬ 
tionally  dependent  on  (ka)m  and  on  the  low-frequency  limit 
of  the  structure  factor  U’(  W),  with  W  =  pvfb/2)  =  w(b/2a)m 
as  the  volume  fraction  of  impenetrable  statistical  particles 
with  diameter  b  2.  2a,  i.e..  in  general,  each  dielectric  particle 
was  visualized  as  having  a  transparent  coating  of  thickness 
(6/2)  —  a.  The  present  paper  applies  the  general  theory'1-4  to 
derive  the  leading  X-dependent  terms  of  t?r  and  tja;  these  de¬ 
pend  explicitly  on  (ka)2  and  a/6  for  all  cases  and  on  appro¬ 
priate  correlation  integrals  A'(w). 

The  correlation  aspects  of  the  distribution  that  we  consider 
are  determined  by  the  statistical-mechanics  radial  distribution 
function5  /(R)  for  impenetrable  particles  and  are  exhibited 
explicitly  as  simple  integrals  over  all  R  of  the  total  correlation 
function  F  =  f  —  1.  The  integrals  for  spherical  and  slab 
particles  are  of  the  form  J  FRndR  (moments  of  F),  but  cy¬ 
lindrical  particles  also  involve  /  FUnRidR.  These  can  all  be 
evaluated  numerically  from  existing  statistical-mechanics 
results  or  approximations-'’  1,1  for  /.  We  obtained  explicit 
closed-form  approximations  before1  17  for  the  integrals  that 
arise  in  the  IT  set  and  also  used  the  required  A  integral  for 
spheres  in  a  related  development111  for  large  kb;  for  slabs,  we 


obtain  both  the  IF  and  ,V  integrals  from  our  earlier  Laplace 
transformation14  of  the  exact  Zernike-Prins  result1-1  for/.  For 
cylinders,  we  may  use  the  virial  expansion  for  /  to  consider 
some  of  the  properties  of  A'. 

In  the  following,  for  brevity,  we  use,  for  example,  form  (4; 
113)  to  indicate  Eq.  (113)  of  Ref.  4,  as  well  as  essentially  the 
same  notation  as  before. 1  4  W’e  generalize  the  earlier  multi¬ 
ple-scattering  electric-dipole  approximations  for  rf2  =  f  given 
collectively  in  form  (1:44)  by  including  scattering  (and  mul¬ 
tipole  coupling)  by  magnetic  dipoles  and  electric  quadrupoles 
(for  spheres  and  cylinders).  For  slabs  and  6  =  2a  (minimum 
separation  of  slab  centers  equal  to  slab  thickness),  the  explicit 
approximation  for  r/2  =  e  reduces  to  f  -*  t'  if  w  —  1,  as  required 
from  physical  considerations:  The  particles  occupy  all  space. 
The  limit  w  -  1  is  not  realizable  for  identical  spheres,  and  we 
take  w  S  u-d  *  0.63,  with  u-d  as  the  densest  random  packing 
introduced  earlier1-17  to  define  the  amorphous  solid.  How¬ 
ever,  our  explicit  approximation  for  e  for  spheres  also  reduces 
to  ('  as  w  reduces  to  1;  we  regard  this  as  consistent  with  the 
approximations  involved  in  scaled-particle8  and  Percus- 
Yevick5-7  statistical-mechanics  theory  and  with  the  closure 
approximation  used  in  the  multiple-scattering  theory.-1-4  For 
cylinders,  we  take  w  S  trd  *  0.84.  as  before.1  Were  an  anal¬ 
ogous  closed  form  available  for  the  .V-integral  for  this  case, 
we  would  expect  the  corresponding  approximation  for  t  to 
show  the  same  behavior  for  the  nonrealizable  limit  w  -*  1. 

The  present  application  of  the  general  theory11-4  to  larger 
kb  than  before1-2  plus  the  recent  applications18  to  large  kb 
provide  simple  forms  that  explicitly  display  the  functional 
dependence  on  all  key  parameters  for  many  practical  appli 
cations.  Thus,  in  these  ranges  of  kb,  elaborate  machine 
computations  are  no  longer  required,  and  the  results  help  to 
delineate  the  fundamental  physical  processes. 

PRELIMINARY  CONSIDERATIONS 

For  a  slab-region  distribution  and  a  normally  incident  wave 
0P-iw/  (representing  either  the  electric  or  the  magnetic 
component)  we  write 
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<j>  =  ee‘kz,  e-z  =  0,  k  =  2?r/A  =  2jtjj,./A,  (1) 

with  T)e  as  the  index  of  the  embedding  medium.  The  corre¬ 
sponding  bulk  coherent  propagation  coefficient 

K  =  ki)blrie  =  ki),  i /'-’  =  «  (2) 

is  to  be  expressed  in  terms  of  p  and  F  for  pair-correlated 
particles  specified  by  their  isolated  scattering  amplitudes  g(r, 
i ).  The  normalization  for  g  is  such  that  for  lossless  parti¬ 
cles 


-Re  e*g(z,  z)  =  -Reg  =  ,W|g(r,  z)|2,  (3) 

with  Jl  as  the  mean  over  all  directions  of  observation  t.  The 
corresponding  known19  scattering  coefficients  an  (which  may 
represent  two  sets)  are  normalized  by  the  form 

g  =  Y.an<  an  =  a„(c\  x),  ('  =  (p/ee,  x  =  ka.  (4) 

In  addition  to  the  dependence  on  t’  (the  relative  dielectric 
parameter)  and  on  x  (the  normalized  radius  or  half-width), 
the  coefficients  depend  on  the  dimensionality  (m)  of  the 
problem  and  on  the  choice  of  field  component  for  m  =  2  (i.e., 
on  whether  the  electric  polarization  is  lateral  or  transverse  to 
the  cylinder’s  axis).  We  obtain  results  for  the  bulk  relative 
dielectric  parameter  e  in  the  form 

t  =  +  (c  +  it,  =  <5 1  +  x2Gc  +  ixm6s,  (5) 

where  the  set  6  is  independent  of  x  =  ka.  The  forms  for  6 1 
and  Gs,  corresponding  to  multiple  scattering  by  electric  di¬ 
poles,  were  discussed  before12  in  detail.  Now  we  obtain 

Sc 

From  Rayleigh's  results  for  spherical  dipoles,20  the  first 
approximation  for  sparse  uncorrelated  distributions  corre¬ 
sponds  to 


Vr  ~  1  =  -cg/2,  (r  =  ij/j2;  c 


i47rp 


i4p 
k 2  ’ 


i2p 

k 


(6) 


In  first  of  the  papers  cited  in  Ref.  1,  for  lossless  small  ('  =  1, 
we  multiplied  Im  ij*  by  the  statistical-mechanics  packing 
factor  W  to  obtain  the  appropriate  r),  for  the  correlated  case. 
The  complete  fi  for  spheres  and  slabs  was  given  by  Maxwell,21 
and  for  cylinders  for  both  polarizations  by  Rayleigh.22  We 
obtained  tj  and  «,  from  forms  (3:74)  and  (4:52) 


T)2  -  1  =  -cG,  rj2  =  (,  (7) 

with  G  as  a  multiple-scattering  amplitude.  This  form  with 
G  =  g  was  obtained  originally  by  Reiche23  and  by  Foldy24  for 
spherical  cases,  and  Lax25  derived  the  form  in  terms  of  a  more 
general  amplitude  than  g.  The  function  G  that  we  require 
is  discussed  in  detail  in  Refs.  3  and  4.  In  particular  for 
spheres,  cylinders,  and  slabs,  respectively,  the  systems  of  al¬ 
gebraic  equations  (4:113),  (3:92),  and  (3:179)  determine  t  = 
tj2  functionally  in  terms  of  an  and  F  for  arbitrary  t'  and  ka  = 
x.  Before,  we  kept  only  the  electric-dipole  coefficient  a  i  in 
these  systems  and  considered  only  the  leading  terms  of  their 
imaginary  and  real  parts,  of  order  xm  and  x2m,  respectively, 
for  lossless  scatterers.  Now  we  include  terms  to  order  xm+2 
for  the  refractive  and  absorptive  effects. 

Thus,  for  spheres  with  x  small,  the  electric  dipole  approx¬ 
imates19 


1563 


ix3(('  —  1)  f  _  x23(2  -  (')’ 
e'  +  2  5(T  +  2) 


xfi2(T  -  l)2 
3(T  +  2)2 


(8) 


where  the  next  terms  are  of  order  ix1  and  x8.  The  corre¬ 
sponding  magnetic  dipole  (a  !m)  and  electric  quadrupole  are, 
respectively. 


Qim 


i’xs(f'  -  1) 
30 


a  2 


ixHf'-  1) 
6(2f'  +  3)  ’ 


(9) 


where  the  next  terms  are  proportional  to  ix7,  ix9,  and  x10. 
The  resulting  amplitude  for  Rayleigh’s  approximation  (6)  is 
g  *  a  i  +  +  a2.  (We  use  an,  a\M  for  the  earlier4  bn, 

Cl) 

For  cylinders19  and  polarization  transverse  to  the  axis,  we 
haveg  ^  ao  +  Qi  +  Q2  with  dominant  dipole 


« l 


iirx2(<'-l)  x2[3  +  <’  —  4(f'  —  1  )L\ 

2(«'+l)  "  8(«'+l) 


X  47T2(f/  l)2 

8(e'+l)2 


(10) 


where  L  =  ln(2/xc')  with  c'  =  1.781  .  . .  The  next  terms  are 
of  order  ix  6  and  x6.  We  also  retain 


a(> 


j  7TX  4  ( f '  -  1) 

32 


a  2 


!JTX4(('  ~  1  ) 

16(e'  +  1) 


(11) 


where  the  next  terms  are  of  order  ix6,  ix8,  and  x8.  For  cylin¬ 
ders  and  polarization  along  the  axis,  we  use  g  ~  alt  +  at, 
with 


a0'- 


ii r2(f'  -  1) 


1  -  —  [3 -e'-(c'-  1)4L] 


x4tt2(('-  U2 
16 


(12) 


where  the  next  terms  are  of  order  ixfi  and  x6.  In  addition,  we 
keep 

d]  *  ijrx4(e'  —  1)/16  (13) 


and  ignore  ix6,  ix8,  and  x8  terms. 

For  normal  incidence  on  slabs,  g  =  a0  +  a  i,  as  discussed  for 
(3:193).  The  dominant  term  is 


a0  =*  ixW  -  1) 


x2(2('  -  1) 
3 


x2(t'  —  1  )2,  (15) 


where  the  next  terms  are  of  order  ixr>  and  x4.  We  also  re¬ 
tain 

a  i  *  ix:1(f'  —  l)/3,  (16) 

but  not  ix  \  ix7,  and  x6  terms. 

DISTRIBUTION  OF  SLABS 

From  form  (3:177)  we  have 

G  =  An  +  A\,  c  =  i‘2p/k  =  iic/x,  w  -  p2a  =  pr(a),  (17) 
where 

A)  =  <i( i(  1  +  AoHo  +  A  i  H  i/»j), 

A |  =  17-0 1  ( 1  +  AoHi/t)  +  Ai’Hn/jf2),  (18) 

with  // 1 1  =  c  +  H 11,  and  ‘Hn  and  H  \  given  as  functions  of  k.  77. 
and  F  in  form  (3:177).  This  algebraic  system  is  valid  ’  for  all 
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ka  =  x,  but  we  consider  only  forms  (15)  and  ( 16)  for  a„  and 
the  corresponding  leading  terms  of  the  correlation  integrals 
Hn  =  //„  +  iS„.  We  have 

/Vo  *  2p  f  *  FdR  +  ik‘2p  C  '  FRdR  =  IF  -  1  +  ixN. 

J  o  Jo 

H  ,  =*  —ixr)N,  (19) 

where  the  next  terms  are  0(.cJ).  Substituting  into  the  ap¬ 
propriate  form  (7)  yields 

e  =  [1  -  ca(i(l  +  a o // o ) ] ( 1  -  cai)  (20) 


Hu  =  2irp  jFRdR  +  (4p  J  F  Mc'kR/2)RdR  -  IF  -  1  +  i.V. 

H „  '  -i(rin2irp/nn)  FrdR  =  -irin(\V  -  l)/nrr, 

(25) 

where  the  next  terms  are  O(x-).  We  obtain  form  (20)  toO(x-), 
from  which 

r  =  1  +  uh  +  x-wh-(  1  +  2  w  +  4M)/8  +  ijrx'2S-u'U’/4. 

=  (26) 

Here 


to  O(.v-).  Thus  in  the  corresponding  form  (5) 

f  =  1  +  in 5  —  jt2tri-’[2  —  tc  +  3NJ/3  +  ixfi-uU’, 

b  =  e‘-l,  (21) 

we  require  only  the  functions  IF  and  .V. 

The  rigorous  pair  function  pf  for  one-dimensional  impen¬ 
etrable  statistical-mechanics  particles  of  width  b  was  given 
by  Zernike  and  Prins.1:1  and  its  Laplace  transform  was  ap¬ 
plied14  in  the  development  of  a  residue  series.  From  the  La¬ 
place  transform  of  f  as  in  Ref.  4,  we  obtain  the  moments  of  F 
=  f  —  1  by  a  Taylor-series  expansion: 

Fu  =  C  '  FdR  =  6(-2  +  W)/2  =  -5F„. 

VV  =  ph  =  pp(b/ 2)  =  ub/2a, 

F i  =  '  FRdR  =  -bM  1  -  4W/3  +  W/21/2  =  -fe-’F,, 

(22) 


IF  =  1  +  2irp  j  FRdR  =  1  +  2xpF\  =  1  —  8 WF j , 

W  =  trp(6/2)-  +  pc  (5/2)  =  w(b!2a)2,  (27) 


and 


M  =  L  —  vS  12 

=  ln(5/a)  +  IF  ln(2/c'kb)  —  27rp  J  F  In{R/b)RdR 
-  ln(6/a)  +  WLh  +  8WT,, 

F,  =  -  f  F(ln  u)udu,  (28) 

with  FiR)  =  F[R/b\  =  F|u]  and  Lh  =  ln(2/c'fe5). 

We  may  evaluate  T,  and  F/  numerically  by  using  tabulated 
values  of  f  or  the  original  integral-equation  approximations 
in  the  computing  routine.  To  first  order  in  W.  we  use  the 
virial  expansion:  F  =  -1  for  u  <  1, 


F 


8W 

ir 


u 

2 


1  <  u  <  2,  (29) 


where  the  notation  is  the  same  as  in  Ref.  18.  (The  identical 
results  follow  on  integrating  the  virial  expansion  of  F  in  powers 
ofW.)  Substituting  into  11"  =  1  +  2pFuand  N  =  (2p/a)F,.  we 
obtain 

IF  =  ( 1  -  W)-.  S  =  -(b/a) W(  1  -  4W/3  +  W-72), 

(23) 

where  IF  was  obtained  earlier1  from  the  rigorous  Tonks 
equation  of  state16  (which  also  follows  from  scaled -particle 
theory6)  by  using  statistical-mechanics  theorems. 

From  elementary  physical  considerations,  if  b  =  2a  (mini¬ 
mum  separation  of  particle  centers  equal  to  particle  width), 
then  for  to  -■*  1  (the  limit  of  a  uniform  slab)  we  require  that  e 
=  ('.  Since  IF  Oanc’  N  --*  -’A,  the  result  (21 )  reduces  tor 
=  as  required. 

More  generally  for  b  =  2a,  the  bracketed  function  in  Eq. 
(21 )  reduces  to 

|  |  =  2-  7ic  +  Sw2  -  :\u  \  5  =  2a.  (21') 

DISTRIBUTION  OF  CYLINDERS 

For  polarization  along  the  axes,  to  the  orders  of  accuracy  in¬ 
dicated  for  forms  (12)  and  (13),  we  follow  the  development  for 
slabs  with 

c  =  i4p/k -  =  w  =  f/n  a2  =  pc(a),  (24) 

where  H\  \  =  (c  +  ’H»  +  H  >)/ 2  in  terms  of  the  correlation  in¬ 
tegrals  in  form  (3:71).  We  have 


and  F  =  0  for  u  >  2.  A  closed-form  approximation  of  IF  (and 
consequently  of  F{ )  was  derived  earlier1  by  differentiating  the 
scaled-particle  equation  of  state.6  Thus  we  obtained1 

VF  =  (1  —  W  )•*/ ( 1  +  W),  (30) 

from  which  IF  =  1  —  4  W  +  7  W-  +  . . . .  The  rigorous  virial 
expansion  to  0(  W-)  is 

IF  =  1  -  4W  +  y/3  12 W2/tt  *  1  -  4W  +  6.6159WL  (30') 

For  the  unrealizable  value  W  =  1,  the  closed-form  IF  vanishes; 
a  comparable  approximation  of  M  would  reduce  to  -:V4  for  b 
=  2a  and  W  =  1  in  order  for  e  to  equal  The  corresponding 
moments  are  then  F,  =  V8andF/  =  —  (In  2)/2. 

For  polarization  transverse  to  the  axes  we  use  G  =  .40  +  .4 1 
+  A  i  with  A„  satisfying  the  system  (3:89).  From  the  solution 
(3:90)  and  the  corresponding  form  of  rj-,  in  terms  of  a„  and  H „ 
of  form  (3:91),  to  the  accuracy  required  for  present  purposes, 
we  work  with 

— («  -  l)/c  =  a0  +  f.4|  +  (-’a s(l  +  .4  ic/2)-’;  (31) 

.4,  =  a,/(l  -  a,  /Vi,). 

2/Vn  =  c  +  (IF  —  1I+i'.V-i<()F-  1)/2jt. 

where  (.4  i  is  a  multiple-scattering  coefficient  that  includes 
all  electric-dipole-dipole  coupling.  The  coefficient  a (>  (es¬ 
sentially  the  magnetic  dipole)  is  uncoupled,  and  the  multiplier 
of  the  electric  quadrupole  a-j  includes  all  orders  of  electric- 
dipole  coupling  to  the  required  accuracy. 

Using  approximations  (10)  and  ( 1 1 ).  we  obtain  initially 
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2u(f'-l)  u'd 

(,  =  1  + - =  1  +  — , 

l  +  1)  I) 


•D  =  1  +  (1  if)&-  &  =  ('—!.  (32) 


with  which 


x2U'd  |3  +  (’  —  d\4M  +  <i ( M’  —  1  )| 


fi(f,  +  l)’l  ,  iirx-S-u'U’  lnnl 
(i  “  }  "4"  ,  *  (33) 

2(c'  +  1)  |  8 »- 

For  b  =  2a  and  the  unrealizable  value  w  =  1,  we  have  '£)  =  1 
and  <i  =  ( and  the  result  for  IF  =  0  and  M  =  —%  is  again  e  = 

For  comparison  with  form  (26),  say,  t/.we  have  to  0(52)  for 
the  present 

(  *  1  +  ud  —  it  (l  —  u)d2/2  +  x-u  d-(  1  +  2it’  +  4 M  +  1F)/16 

+  nrx-S-w  UV8.  (34) 

The  corresponding  bipolarization-  is 

(i  —  e,  «  u(l  -  u  )6-/2  +  x'-ici-ll  +  2ic  +  4M  —  U’)/16 

+  nrx-5-w'U’/8,  (35) 

and  we  may  obtain  higher-order  terms  in  <5  from  Eq.  (33). 


DISTRIBUTION  OF  SPHERES 

For  spheres  we  use  form  (7)  with 

G  =  A i  +  Aim  +  t42.  c  ~  i4irp/k2  =  t3u7x3, 

ic  =  p4rra:V3  =  pv(a),  (36) 


.  .  3u>U'  -  1)  ,  ,  wd 

(1  =  1  + - —  1  +  ’ 

2  +  f'  —  w(i  —  1 )  » 


»  =  i  - 


(1  —  u>)6 


d-f'-l,  (40) 


with  which 


2  A  I  1  f2  ~  N  to 
<  =  f i  —  x2wd  — -  - H —  2  +  — 

I®2]  5  9  l  5 


(2c,  +  3)211  ix3262wW 


1+— - -  + - - — •  (41) 

10  l  5(2e'  +  3)Jj  9S2 

The  magnetic-dipole  contribution  x‘2wd(i  atM(i  shows  that 
the  effects  of  the  function  in  brackets  in  Eq.  (37)  have  can¬ 
celed. 

By  differentiating  the  scaled-particle  approximation6  for 
the  equation  of  state,  we  showed  that1,17 

1F  =  (1  -  W)V(1  +  2W)2,  (42) 

which  also  follows  from  the  Percus-Yevick  approximation.7 
The  first  moment  F i  obtained8  from  the  Wertheim-Thiele 
solution  of  the  Percus-Yevick  integral  equation7  gives  the 
closed  form 


b  1  +2W 


_  W  W* 

5  10,  ' 


Although  the  physically  realizable  range  corresponds  to  W  < 
Wd  ~  0.63,  we  see  that,  for  b  -2a  and  the  unrealizable  value 
tr  =  1,  it  follows  that  U’  =  0  and  N  =  9/5;  then  c  =  e ',  as  was 
discussed  for  slabs  and  cylinders. 

For  comparison  with  forms  (21),  (26).  and  (34),  we  have  to 
0(^2)  for  the  present  case  of  the  sphere 

c  *  1  +  u’d  -  «d  -  u  )i 52/3  +  x  -u  S2(6  +  3 «■  -  5JV)1 1/(15)-’ 

+  ijc32(5VU>/9.  (44) 


with  the  A's  satisfying  the  system  (4:113)  in  terms  of  the  iso¬ 
lated  coefficients  an  and  the  correlation  integrals  H  „  of  form 
(4:80)  or  (3:148).  Introducing  the  low-frequency  forms  (3: 
149),  we  solve  the  system  and  obtain 

-(f)2  -  l)/r  =  V2A 1  +  f)2aiA/[l  +  A\ct)I(t)  +  l)]2 
+  f)4a2(l  +  A  ic3/5)2; 

A  i  =  a,/(l  -  a\‘H\\), 

3 Hu  =  2 c  +  2 -H0  +  H2  *  2c  +  2(11’  -  1) 

+  i/V(2+ r)2/5)/x,  (37) 


For  b  =  2a,  the  function  in  brackets  reduces  to 

|  )  =  3(2 -5ir  +  4u-2-u  :')/(1  +  2u  ).  (44') 

BULK  INDEX  OF  REFRACTION 

We  write  forms  (21),  (26),  (33),  and  (41)  collectively  as 

«  =  t]  +  wd2P(x2 )  +  iwd-S(xm),  &  =  c'  —  1,  (45) 

wh  re  P  and  S,  proportioned  to  x2  and  xm,  are  obtained  by 
i  ction.  The  k -independent  term 


with 

IF  =  1  +  47rp  £  FR2dR  =  1  +  4ir pF2  =  1  -  24WF2, 

W  =  p4fr(6/2)  V3  =  pv(bf 2)  =  u  (b/2a  )-'*,  (38) 

and 

N  =  -4irpa  FRdR  =  -4npaF]  =  24(a/h)WF|.  (39) 


Here  rj2.A \/a\  includes  all  electric -dipole-dipole  coupling,  and 
the  multipliers  of  the  magnetic  dipole  a, m  and  the  electric 
quadrupole  a2  incorporate  multipole  coupling  with  all  orders 
of  electric  dipoles  to  the  required  accuracy. 

(/sing  approximations  (12)  and  (13).  we  obtain  initially 


Ci  =  1  +  !c<5/(1  +  dD),  D  =  (1  —  u)Q, 

Q ,  =  Q21  =  0,  Q2<  =  '/2,  Q:<  =  ‘/:,  (46) 

represents  special  cases  of  the  result  for  ellipsoids.2  The 
corresponding  S  for  ellipsoids  is  also  known,2  and  form  (45) 
holds  for  all  small  dielectric  particles  (discounting  a  resonant 
multipole). 

The  corresponding  bulk  index  of  refraction  may  be 
written 


f)  =  |c ,  +  wdHP  +  /.S)| l/-  *  f),  +  wdMP  +  /S)/2f?i,  (47) 
with  f),  =  y  c i-  More  generally,  we  write  t  -  fr  +  ii,  and  t)  = 

f?r  + 


1/2 


M  =  [  f  r  2  +  (48) 
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In  terms  of  A  =  c  —  1,  we  obtain  form  (2:3),  so  that  part  of  the 
earlier  development  is  appropriate. 

In  particular,  if  we  retain  c ,  to  0(5:i)  and  P  and  .S’  to  0(5"), 
then 

U  -  1  )/w  =  A/ir  =  5(1  -  5D  +  &-D-)  +  52<P  +  iS)  (49) 

is  correct  to  0(53),  0U-5-),  and  i'0(xm52).  For  complex  5  =  5r 
+  (5,  vve  construct  A  =  Ar  +  i  A,  essentially  as  for  form  (2:25) 
with  the  earlier  iS  replaced  by  iS  +  P  but  retain  only  the 
leading  term  of  the  earlier  5:i  contribution.  Thus 

Ar/u  =  5r(l  -  hrD  +  5r2D2)  +  5,2D<  1  -  3 5rD) 

-  2 brhjS  +  (5r2  —  5 ,  -)P  (50) 

and 

X/w  =  5,|1  -  2 brD  +  <35r2  ~  5, -)D-\ 

+  (5r2  -  o, -).S  +  2 5r5,P.  (51) 

The  terms  in  P  account  for  the  O(Ar-)  corrections  indicated  for 
forms  (2.22)  ff.  If  5,  =  0,  then  Ar  is  independent  of  S,  and  A, 
of  P. 

Similarly  in  terms  of  «<  =  tj'  —  1  with  tj'  =  7)p/7j,, ,  we  write  5 
=  r/'-  -  1  =  e(2  +  id  and  express  ?/  -  1  to  O(r'),  0 (ar and 
(0U "•()-*)  as 

(17  -  l)/tc  =  i-  +  i'-’(. 4  +  4P  +  i'4S’)/2  +  i>:1B/2, 

.4  =  1—  («•  +  4/1),  H  =  -(1  -«)(«+  4ZJ)  +  SD~.  (52) 

For  complex  r  =  vr  +  it’,  we  construct  r/  =  rjr  +  it),  essentially 
as  for  form  (2:28).  Thus  the  refractive  contrast  rjr  —  1  is  given 
by 

U7 r  -  ll/«  =  Ml  +  iy.4/2  +  2)  -  t'.HA  +  3i'rB)/2 

-  4i’,rrS'  +  2(i’r-  —  i’,-)P.  (53) 

and  the  net  attenuation  is  determined  by 

tj./k  =  i', [1  +  iy.4  +  ( 3 i’r -  -  2) 

+  2dv2  ~  i',-).S  +  4 1', i’rP.  (54) 
If  the  .-catterers  are  lossless  (/',  =  0),  then  r\r  does  not  depend 
on  N.  and  77,  is  independent  of  P.  The  present  forms  (50)— (54) 
hold  for  all  small-spaced  particles  to  the  indicated  accuracy. 
Generalization  of  b  and  S  to  aligned  ellipsoidal  particles  and 
cocentered  nonsimilar  ellipsoidal  exclusion  regions  are  given 
in  Ref.  2. 

For  measurements  in  which  the  parameters  r0  and  7j()of  the 
embedding  medium  are  varied,  we  normalize  (j,  and  jjt,  with 
respect  to  ep  and  tjp  instead  of  and  77,,.  See  Ref.  2  for  de¬ 
tails. 

The  values  of  I).  P,  S.  A,  and  B  for  forms  (49)— (54)  are  given 
by  the  following  for  the  special  cases  at  hand.  For  slabs, 

I)  i=0,  /4i  =  1  —  ic,  B|  =  — ir(l  —  «’),  (55) 

P,  =  -x2(2  -  W  +  3JV)/3,  .S,  =  j  IP,  (56) 

with  IP  and  ,V  as  in  forms  (19)  and  (23).  For  cylinders  and 
lateral  polarization,  1 1)2i,  A->i,  B-ti  1  equals  |D,,  A 1,  B|],  and 

P-,i  =  x  -’( I  +  2  u-  +  4M)/8,  S2i  =  jr-jr  ll’/4.  (57) 

For  cylinders  and  transverse  polarization, 

I)-,,  =  (I  -  11 9/2.  A  ,,  =  -(1  -  a  )  =  -A,,. 

B_,,  =  — /c(  1  -  a  )  =  H  >i,  (58) 


P-2,  =  x2(l  +  2  w  +  4  M  +  1P)/16  =  P-ut2  +  x-lP/16, 

S-n  =  ttx  2‘lP/8  =  S21/2,  (59) 

with  M  and  ‘IP  as  in  forms  (25)-(30).  For  spheres, 

U3=U-w)/3,  =  — (1  —  tc)/3. 

Pa  =  -(1  —  ic)(4  +  5u.')/9,  (60) 

Pa  =  x2(6  +  3w  —  5/V)ll/(15)2,  S:<  =  2x"W/9,  (61) 

with  N  and  VP  as  in  forms  (38),  (39),  (42),  and  (43). 

For  cylinders,  the  values  of  t'  and  17'  for  the  lateral  and 
transverse  cases  may  differ,  and  the  corresponding  birefren- 
gence  (7)7,  -  t)(I  )  will  then  display  intrinsic  as  well  as  form  ef¬ 
fects.  See  Ref.  2  for  details.  The  relations  A 21  =  —A 21  and 
P21  =  B21  simplify  considerations.  In  particular,  if  there  is 
a  common  v ,  then 

77/  —  7),  =  ui>2\Ai  +  Pi  -  x2‘)P/16  +  iS/) 

=  wv2[l  -  w  +  x 2(  1  +  2 u>  +  M  -  )P)/16 
+  i  7rx  2‘VP/4  ]  =  i’HR  +  il).  (62) 

Then  the  birefrengence  corresponds  to 

Re( 77/  -  77,)  =  (17  2  -  »i2\R  -  2 tv  15 /  (63) 

and  the  dichroism  to 

Im(77/  -  7), )  =  2 vrv,R  +  (17  2  -  c,2)/.  (64) 

For  experiments  in  which  ije  is  varied,  we  introduce  the  vari¬ 
able  £  =  (7 ),.  -  Vpr)/Vpr  and  the  constant  n  =  r)pJr\pr  to  con¬ 
struct  v  -  (-£  +  ip)/(  1  +  £).  For  small  £.  we  have 

Re(7)h/-  Vbi)/ripr*  [(£2-m2>R  +  2£h/|(1  -  £),  (65) 

Im(77W-  77s,)/v  *  |-2£^P  +  (£2-t72)/|(1  -£),  (66) 
which  generalize  the  result  in  form  (1:67). 
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